Growing correlations and aging of an elastic line in a random potential 
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We study the thermally assisted relaxation of a directed elastic line in a two dimensional quenched 
random potential by solving numerically the Edwards- Wilkinson equation and the Monte Carlo 
dynamics of a solid-on-solid lattice model. We show that the aging dynamics is governed by a 
growing correlation length displaying two regimes: an initial thermally dominated power-law growth 
which crosses over, at a static temperature-dependent correlation length Lt ~ T 3 , to a logarithmic 
growth consistent with an algebraic growth of barriers. We present a scaling arguments to deal 
with the crossover-induced geometrical and dynamical effects. This analysis allows to explain why 
the results of most numerical studies so far have been described with effective power-laws and also 
permits to determine the observed anomalous temperature-dependence of the characteristic growth 
exponents. We argue that a similar mechanism should be at work in other disordered systems. We 
generalize the Family- Vicsek stationary scaling law to describe the roughness by incorporating the 
waiting-time dependence or age of the initial configuration. The analysis of the two-time linear 
response and correlation functions shows that a well-defined effective temperature exists in the 
power-law regime. Finally, we discuss the relevance of our results for the slow dynamics of vortex 
glasses in High-Tc superconductors. 

PACS numbers: 74.25.Qt,75.60.Ch,64.70.kj,64.60.Ht 



I. INTRODUCTION 

Disordered elastic manifolds play an important role in 
a variety of physical systems. Growing surfaces inter- 
faces in domain growth phenomena and phase separa- 
tion^ and cracks in brittle materials^ are usually viewed 
as elastic objects. Interesting one dimensional realiza- 
tions are polymers 4 and vortex flux lines in type II su- 
perconductors.™ While polymers can wrap, vortex flux 
lines are preferentially directed along the applied mag- 
netic field. 

A number of systems involving one dimensional elastic 
manifolds display glassy features. An ensemble of inter- 
acting polymers forms a polymer melt that undergoes a 
glass transition^ The competition between vortex repul- 
sion and their pinning to randomly located impurities 
also leads to glassy phases in superconductors.— While 
in the former case disorder is self-induced, in the latter 
the effect of the impurities is mimicked by a quenched 
random potential. 

The relaxation dynamics of a model of layered high 
T c superconductors was recently studied in Ref. 0. The 
magnetic vortices are considered to be directed along one 
direction and they can move in two transverse directions. 
Conventionally, this is a (d = 1) + (N = 2) dimensional 
model. The dynamics following a rapid quench from the 
liquid state to a set of control parameters in which the 
equilibrium state is expected to be the so-called 'vortex 
glass— was monitored. By using numerical simulations a 
slow out of equilibrium relaxation typical of glasses was 
found. In this system, diffusive aging of the averaged 



two-time roughness, (w 2 ), and displacement, (B), was 
observed. The aging was characterized by a multiplica- 
tive scaling form £ 2 ^ (t w ) f w 2 ^ B (£(t) / £(t w )) , where £(t) is 
a growing characteristic length. The scaling functions 
f w 2 and /b are different but the exponent £ does not 
depend on the observable (w 2 or B). Besides, the sys- 
tem of interacting lines was perturbed in order to com- 
pute the integrated linear response, and a diffusive aging 
was also found, characterized by a scaling function of 
the type £ 2< > (t w ) f x (£(t) / £{t w )) with the same exponent £ 
and growing length £(t) found for the unperturbed two- 
time observables (roughness and displacement). It was 
also shown that correlation and response functions can 
be related by a modified fluctuation-dissipation relation 
after removing the 'diffusive' contribution, i.e. the fac- 
tors £ 2 ^{t w ). This violation of the fluctuation-dissipation 
relation can be characterized by a finite effective tem- 
perature, Toff- 7,8 Furthermore, within the time-window 
that was numerically explored, the growing length £(t) is 
well-described by a power-law, £(t) ~ i 1 / 2 . This could 
well be a pre- asymptotic regime after which the growing 
length crosses over to the expected activated dynamics 
logarithmic growth (when disorder free-energy barriers 
scale as A(L) ~ lr with ip > 0). It is worth to mention 
that although aging in high 9 and low 10 T c superconduc- 
tor samples has been reported in the literature, a detailed 
comparison to the results listed above remains to be done. 

The model system used to describe vortex lines in high 
temperature superconductors contains several contribu- 
tions coming from different energy scales. The main con- 
tributions are: the elastic energy of the individual lines, 
non-linear terms in the elastic energy, the interaction be- 
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tween the lines, and the quenched random pinning poten- 
tial. It is clearly important to establish whether all the 
above listed contributions to the energy are necessary 
to get such aging behavior or whether similar features 
arise when some of the terms - non-linear contributions 
to the elasticity, interactions between the lines, random 
potential - are switched off. We review below the out of 
equilibrium dynamics of related models with and without 
this type of interactions. 

The first studies of the out of equilibrium dynamics 
of directed elastic manifolds in a quenched random envi- 
ronment focused on mean-field models in which the trans- 
verse space is infinitely dimensional, N — > oo, each trans- 
verse coordinate has infinite length, M — > oo, and the 
manifold has finite dimension, d, and infinite length in 
all directions, L — > oop^ 1 - This model includes only two of 
the energetic contributions listed above viz. elasticity and 
quenched disorder and neglects non-linear terms and in- 
teractions between the elastic objects. A dynamic phase 
transition separating a liquid (high-temperature) phase 
from a glassy (low-temperature) one was found for all in- 
ternal dimension d including d — 0. 11 Different aging dy- 
namics characterize the low-temperature phase depend- 
ing on the short or long range character of the random 
potential correlations. In both cases the aging regime 
lasts for ever (after having taken the N — > oo, M — ► oo 
and L — > oo limits): there is no multiplicative factor 
^{t) in the averaged two-time observables. For short- 
range correlated potentials the displacement, roughness 
and linear responses scale as f w 2 b tX (t(t) / £(t w )) and a 
single finite effective temperature exists. 

Soon after Barrat 1 ^ and Yoshino 13 ' 14 studied a solid- 
on-solid model of a single directed elastic line relaxing 
with Monte Carlo dynamics on a disordered substrate 
with one (N = 1) transverse direction. This system mod- 
els elasticity in a rather extreme way, using a hard con- 
strain, and includes quenched randomness. There are no 
interaction between different lines in this model. These 
authors focused on the very long length limit in an ef- 
fectively infinite box (M L ^> a with a the lattice 
spacing) and found that this lattice model has a similar 
averaged dynamics to the one found later in the vortex 
glass^ below a crossover temperature. The relaxation is 
slow, and the global displacement and linear response 
show non-trivial diffusive aging with multiplicative scal- 
ing. The characteristic length scale also appears to be a 
power-law of time within the numerical time window and 
the exponents z and a = C,/z are temperature and disor- 
der strength dependent with the same qualitative trend 
as in the fully interacting case. More recently, we focused 
on the analysis of the averaged and fluctuating two-time 
roughness of such solid-on-solid lines with finite length, 
L < ooJ^ On the one hand, we found that the aging 
regime stops and crosses over to saturation of the two- 
time roughness and free diffusion of the displacement at 
a characteristic value of the time delay, At = t — t w = t x , 
that grows with the length of the line and smoothly de- 
pends on other parameters in the model. The saturation 



of the two-time roughness is well described by a gener- 
alization of the equilibrium Family- Vicsek scaling^ On 
the other hand, the two-time roughness fluctuations are 
highly non-trivial but can be characterized with a rela- 
tively simple argument by a scaling function. 

The numerical solution to the Langevin equation for 
free lines in N = 2 transverse dimensions 6 and, espe- 
cially, the full analytic solution to the Langevin dynamics 
of finite Edwards- Wilkinson (EW) elastic lines^ in one 
transverse dimensio n 18 i 19 (TV = 1) had also been con- 
sidered, without taking into account interactions, non- 
linear terms and quenched randomness. The results 
suggested that not even disorder is necessary to ob- 
tain similar averaged and fluctuating aging dynamics of 
fully relaxing quantities. The aging scaling and sat- 
uration phenomenon of the averaged two-time rough- 
ness follow the same scaling laws as above with growing 
length £(t) ~ i 1 / 2 , and temperature- independent expo- 
nents z — 2, £ = 1/2 and a — C/ z — 1/4- However, 
the noise-averaged linear response is stationary in this 
'quadratic' model. When it comes to analyzing fluctu- 
ations other differences appear. The distribution of the 
two-time roughness satisfies a similar scaling law as in 
the disordered problem although with a rather different 
scaling function and the linear response simply does not 
fluctuate. 

The effect of non-linear terms has been considered by 
Bustingon y 20 ! 21 who analyzed the relaxation dynamics 
of 'clean' finite length [L < oo) Kardar-Parisi-Zhang 
(KPZ) finest in N = 1 transverse dimensions. In this 
work correlations were analyzed in detail and undergo 
diffusive aging with the KPZ exponents £, a and z. The 
two-time and length scaling of averaged quantities and 
the scaling form of the probability distribution functions 
proposed in Ref. were thus confirmed. 

Finally, once the elastic lines are allowed to wrap, bi- 
ologically motivated dynamic problems can also be ad- 
dressed. A recent study concentrates on the effects of 
confinement on the out of equilibrium relaxation of a 
single polymer chain in two dimensions* 2 ^ a problem of 
relevance for cellular modelling. Another study analyzes 
numerically the effect of randomly applied forces on an 
ensemble of interacting polymer lines, focusing on out of 
equilibrium properties of active matter* 2 ^ There is, cer- 
tainly, much room for further investigations in the bio- 
logical context. 

In this paper we return to the problem of the relax- 
ation of directed elastic lines in the presence of quenched 
randomness. Interactions among different lines and non- 
linear terms are not considered here. Our aim is to com- 
plete the analysis of the averaged two-time observables. 
In order to compare with previous results and to extract 
the universal behaviour we treat two models in paral- 
lel: the Monte Carlo dynamics of the disordered solid- 
on-solid mode l 12 i 13 ' 14 and the Langevin dynamics of the 
disordered Edwards- Wilkinson equation in (1 + 1) dimen- 
sions. We present a careful study of finite-length effects 
on the averaged dynamics of correlation and linear re- 
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sponse functions. We pay special attention to the behav- 
ior of the crossover length between thermal and disorder 
dominated regimes, and discuss finite-size and finite-time 
effects. Moreover, we show by using scaling arguments 
how this length scale determines the relaxation properties 
of the system. 

The organization of the paper is the following. In 
Sect. [U we define the models and we describe the nu- 
merical method and two-time observables on which we 
focus. In Sect. Iffll we measure the growth, roughness and 
dynamic exponents by analyzing the two-time structure 
factor and two-time roughness. We also study in detail 
the growing length £(t) and we confirm the existence of a 
crossover from power-law to logarithmic growth that is, 
however, only seen for sufficiently long lines. Section ITVl 
summarizes the two-time behavior of the averaged rough- 
ness and associated linear response using different initial 
conditions. In Sect. [V] we present scaling arguments to 
explain many of the features observed numerically. Fi- 
nally, in Sect. ["VTl we present our conclusions. 



II. 



THE MODELS 



In this section we introduce the models and the nu- 
merical methods used to study their relaxation. 



A. Lattice model 

A discrete model of a one-dimensional directed elas- 
tic object represents the line as a lattice string of length 
L directed along the y direction i 12 i 13 i 14 i 15 The line can 
move transversely along the x direction on a rectangular 
square lattice of transverse size M = 10 4 3> L 2 / 3, ensur- 
ing the existence of many nearly equivalent, quasi-ground 
states. The line segments, x{y) (y = 1, • • • , L), obey the 
restricted solid-on-solid (SOS) rule \x(y)—x(y—l)\ = 0, 1. 
A quenched random potential V taking independent val- 
ues on each lattice site is drawn from a uniform distri- 
bution in [—1,1]. We use 10 5 — 10 7 realizations of the 
quenched randomness depending on the value of L. At 
each microscopic time step we attempt a move of a ran- 
domly chosen segment to one of its neighbors restricted 
by the SOS condition and we accept it with the heat-bath 
rule. One Monte Carlo (MS) step is defined as L update 
attempts. We use angular brackets to indicate the aver- 
age over thermal noise realizations. We choose two types 
of initial conditions: equilibrium at high temperature ob- 
tained after evolving a random initial condition during a 
sufficiently long time interval at high temperature, and 
equilibrium at zero temperature. 

The lattice model has no finite elastic energy. Elastic- 
ity is modelled in an extreme way; in the absence of dis- 
order all configurations have equal vanishing energy but 
the displacement between neighboring bonds is bounded 
to —1, 0, 1 lattice spacings. 



The control parameters are te mperature, T, and the 
disorder strength, Vo = ^/[V 2 ]y. Here and in what fol- 
lows we use square brackets, [. . . ]y, to indicate an av- 
erage over quenched randomness. The Monte Carlo rule 
implies that the dynamics depend only on the ratio be- 
tween these parameters, Vq/T. In the simulations shown 
here we fixed Vo = 1/3. In the following we use adimen- 
sional time, space and energy scales. 



B. Continuous model 

The SOS model is easy to simulate but it is not sim- 
ply recovered as a limit of better-known continuous prob- 
lems like the KPZ equation 20 ^ or the vortex line model 6 
in which elasticity is modelled in a more realistic way. 
For this reason we also study the disordered Edwards- 
Wilkinson (EW) linei 

The disordered EW equation^ 7 - for a scalar field x rep- 
resenting the height of a surface over a one dimensional 
substrate parametrized by the coordinate y (a one dimen- 
sional directed interface) is 

jd t x{y, t) = cd 2 x(y, t) + F p [x(y, t)} + h[x(y, t)] 

H(y,t), (i) 

where £ is a Gaussian thermal noise with t)) = 
and 



(Z(y,tW,t'))=27T5(y-y')5(t-t') 



(2) 



The parameter c is the elastic constant, 7 is the friction 
coefficient, T the temperature of the thermal bath (in 
energy units) and (• • • ) the average over the white noise 
£, i.e. the thermal average. The Boltzmann constant has 
been set to one, fcg = 1- The term h represents the effect 
of a perturbing field that couples linearly and locally to 
the height, —h(y,t)x(y,t). One can also consider other 
types of perturbation that couple to more complicated 
functions of the height as defined in Sect. Ill Dl The ran- 
dom pinning force F p [x(y, t)] = —d x V{x,y) represents 
the effect of a random-bond disorder described by the 
potential V(x, y), whose sample to sample fluctuations 
are given by 



[V(x,y)-V(x',y')Y 



6(y - y') R(x - x 1 ) , (3) 



where R(u) stands for a short ranged correlator along the 
x coordinate with range 77. The continuous random po- 
tential is modeled by a cubic spline passing through M 
regularly spaced uncorrelated Gaussian number points^ 
We adimensionalize Eq. {T]) by using r j as the unit of dis- 
tance in the x direction, 77*/ as the unit of time, and Vo 
as unit of energy/temperature. The unit of distance in 
the longitudinal direction y can be conveniently taken as 
the layer spacing s of the numerically discretized lapla- 
cian (such a choice is natural when modeling a layered 
material such as a High-T c superconductor with an ex- 
ternal magnetic field applied perpendicular to the oxide 
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planes). In these units, the friction and elastic coeffi- 
cients are equal to one, leaving only two independent 
parameters in the model: the adimensionalized elastic 
constant v — crf/s 2 and the adimensionalized tempera- 
ture T/Vq (from now on, we use T for T/Vq)- Finally, by 
choosing v = 1 we focus exclusively in the temperature 
dependence for a fixed disorder. 

We use a finite-difference algorithm to integrate the 
partial differential equation in which the first and sec- 
ond order partial derivatives are discretized in the usual 
way. (See Ref. for more details on the numeri- 
cal technique.) We typically simulate lines with length 
L = 64, 256, 1024. The time step is t = 0.01. We use 
10 3 noise realizations to evaluate the two-times averaged 
correlations and responses. 

Note that in the continuous model each non strictly 
flat configuration has a finite elastic energy while in the 
lattice model all configurations have the same vanishing 
elastic energy. 



C. Quenched randomness 

The kind of quenched disorder used in both models is 
of the 'random bond' type in the sense that, for a domain 
wall described by our elastic interface model, it does not 
break the symmetry properties of the corresponding or- 
der parameter^ This implies that the two dimensional 
disorder potential locally couples to the interface posi- 
tion and that R(u) in Eq. (J3j> saturates to a constant 
for large distances. The way we generate the disorder 
corresponds to a short-ranged function R(u). The effect 
of long-range correlated randomness has been studied in 
mean-field elastic manifold models^ but we do not con- 
template it here. 



of the lines, a quantity that has been widely used in 
the study of interface dynamics,— but generalized here 
to include the t and t w dependence. The two-time 
roughness is given by 



(w 2 ){t,t v 



Sx(y,t w )Y 



(4) 



where 5x(y,t) — x(y,t) — x(t) accounts for the displace- 
ment of the y-th line segment relative to the center of 
mass, x(t) = L- 1 x(y, t). 

The two-times structure factor is defined 

{S n )(t,t w ) = L {[\c n {t) - c n {t w )\ 2 ] v ) , 

Lc n (t) = £[ a r(y,i)-3f(t)]e- < «»"' (5) 



where q n = 2nn/L with n integer. 

Further insight into the dynamics of out of equilibrium 
systems is given by the linear response function. The 
latter is defined by applying a random time-independent 
force at a time t w on a replica of the system, and by 
computing how this one departs from an unperturbed 
one evolving with the same thermal noise. Concretely, 
the energy contribution of a field conjugated to the dis- 
placement with respect to its mean value is 



^[x(y,i)-Ic(t)] S (y)0(Ai) 



(6) 



s(y) — ±1 with equal probability, ((s(y)}) — and 
((s(y)s(y'))) — 5 y ,y^ with ((. . .)) denoting the average 
over the perturbing field distribution, h is the intensity 
of the perturbation. The associated linear response func- 
tion is 



D. Observables 



f 



(X)(t,t w ) = ([[Sx h (y,t) - Sx(y,t)]s(y)] v ) . (7) 



We study the relaxation after two types of rapid 
changes in the control parameters. The first protocol, 
and the more usual one, describes a quench from high 
to low temperatures. After equilibration at high T the 
system is quenched to a low value of T. The second pro- 
tocol, less usual, consists in first equilibrating the sample 
at zero temperature (i.e. starting in its ground-state), 
To = 0, and then heating the sample to a higher tem- 
perature. In both cases, the line is allowed to relax from 
the quench occuring at time t = until a waiting-time 
t w , when the quantities of interest are recorded and later 
compared to their values at subsequents times t > t w . 

The aging dynamics of elastic lines has been 
initially studied in terms of the averaged two- 
time mean-squared-displacement (B)(t,t w ) = 

The behavior of this 
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x{y,t) - x{y,t w )Y ^ ^ 

quantity is somehow obscured by the motion of the 
center of mass. We thus prefer to focus on the roughness 



Henceforth (. . .) indicates the average over the thermal 
noise and the s(y) distribution. 

In equilibrium the averaged linear response is related to 
the averaged spontaneous fluctuations of the correspond- 
ing observable by the model-independent fluctuation- 
dissipation theorem (FDT) which states 



(w 2 )(At) = 2T( x )(At) 



(8) 



(the Boltzmann constant has been set to one, ks = 1), 
where the At argument implies stationary dynamics. In a 
system relaxing out of equilibrium this relation does not 
necessarily hold. In a number of glassy systems one can 
define an effective temperature^ from the modification of 
the above relation. In the aging regime of elastic lines in 
disorder media the FDT is violated and one constructs 
the modified FDT&&1 3 . 



(w 2 )(t,t w ) = 2T cB (t,t w ) (*)(*,*» 



(9) 
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The two-time dependence of the effective temperature 
is here kept general, ft turns out that in models with 
multiplicative scaling, as the one discussed here, once the 
factors i^(t) have been taken into account, the re-defined 
effective temperature appro aches a constant value (see 
Sect. HVC] and Refs. iliallMIIllil- 



E. Dynamic crossover 

In (I + I) dimensional models of elastic mani- 
folds the glassy phenomenon appears as a dynamic 



crossover 



12,13.14,15.18,19,20 



For all observation times, t bs, that are longer than a 
size, disorder-strength and temperature dependent equi- 
libration time, t eq , the system reaches equilibrium. In- 
stead, for t bs < t eq the relaxation is non-stationary and 
thus occurs out of equilibrium as demonstrated by two- 
time correlations and linear responses that age and by a 
non-trivial relation between them. 

The equilibration time t eq increases by increasing the 
size L of the line, by decreasing the temperature T, and 
by increasing the quenched disorder strength Vq. At fixed 
disorder strength and line length the dynamic crossover 
reflects in a temperature crossover that resembles a phase 
transition at a characteristic temperature T co . At high 
temperature, T > T co , the equilibration time is short, 
disorder is irrelevant and the line behaves as the clean 
EW one at high temperature: the dynamics of a finite 
length line reaches equilibrium meaning that stationary 
dynamics and the FDT hold, and the exponents (de- 
fined and discussed below) are the ones of the clean EW 
line, a = 1/4, £ = 1/2 and z = 2. At lower temper- 
atures, T < T co , the equilibration time is longer than 
the observation time and the dynamics remains non- 
stationar y 12 i 13 i 14 i 15 Quenched disorder modifies the val- 
ues of the exponents and induces a geometrical crossover 
at a temperature dependent length-scale Lt, from a short 
length-scale roughness described by the 'thermal' expo- 
nent (t, to a large length-scale roughness described by 
a 'disorder' exponent (d- We shall discuss the values of 
the exponents and the out of equilibrium dynamics in the 
following sections. 



III. 



GROWTH AND SATURATION 



In this Section we recall the main tools used to analyze 
the evolution of the conformational properties of elastic 
manifolds and we generalize them to take into account 
the effect of the waiting-time. This leads us to present 
the temperature dependence of the crossover length scale 
Lj 1 . 



A. Roughness 

1. Comparison to the initial condition 

Traditionally, the dynamics of elastic manifolds has 
been classified in universality classes according to the be- 
havior of the two-time roughness ([4]) evaluated at t w = 0, 
the initial time, and using a flat initial configuration, 
x(y, 0) = xq*^ The waiting-time being identical to zero 
all two-time observables depend on At = t when com- 
pared to the initial condition. Initially the roughness in- 
creases as a function of At, and at a characteristic time 
t x {L) reaches saturation at an L-dependent value, (w 2 ^). 
This behavior is encoded in the Family- Vicsek scaling 16 
that, in full generality, can be expressed as 



(w 2 )(At) 
lim (w 2 )(At) 



<7(At) , 

f(L), 

h(L), 



(10) 

(11) 
(12) 



with h,g,f three monotonic functions. Consistency at 
At = t x requires h = f o g. In the usually discussed 
space-time scale-invariant cases in which all functions are 
power-laws one has 



(w 2 )(At) 

(wL) 



At 2a 
L z , 
L 2 < 



(13) 
(14) 
(15) 



with a the growth exponent, z the dynamic exponent, 
and £ the roughness exponent. Consistency implies that 
the three exponents are related by 

za = (. (16) 

The values of the exponents are well-known in a number 
of cases. For the EW elastic line the exponents can be 
analytically computed and one finds that a = (2 — d)/A, 
z = 2 and C = (2-d)/2 for d < 2. For the non-linear KPZ 
elastic line the exponents are known only numerically 
for general d. In (1 + 1) dimensions, however, they can 
be computed analytically and a = 1/3, z = 3/2 and 
C = 1/2. 

In the presence of quenched disorder the dependence 
of the asymptotic roughness {w 2 ^) with the length of the 
line undergoes a crossover. For lines that are shorter than 
a temperature and disorder strength dependent value Lt 
the behavior is controlled by thermal fluctuations and the 
relation (fTS")) holds with C = Ct the thermal roughness 
exponent. This exponent is the one corresponding to 
the EW equation, and thus (t = (2 — d)/2 in general, 
and = 1/2 in our (1 + 1) dimensional case. In this 
thermally dominated scale, the dynamics is expected to 
be 'normal' in the sense that lengths and times should 
be thus related by power-laws of the type (fl3|) - (fl5|) with 
the exponents linked by Eq. fTB]). 

For lines longer than Lt, the roughness is dominated 
by quenched disorder and one has that (TT5|) still holds 



6 



though with a different value of the roughness exponent, 
C = Cd- The disorder dominated roughness exponent C,o 
is the one characterizing the geometry of the ground state 
configurations, and it is expected to satisfy Co > Ct in 
general^ In (1 + I) dimensions one has Cd = 2/3. 

Once quenched disorder is present, the time evolution 
is expected to be driven by activation over free-energy 
barriers. If they scale as U(L) ~ L^, the Arrhenius 
law leads to a logarithmic relation between equilibrated 
lengths and times that implies t(£) ~ e T (^/ L r) i ' /T w m 1 
T and Lr some characteristic energy and length scale, 
respectively. One might then expect 

(w 2 )(At)~ln ai "At, (17) 

which for consistency implies 

i> a la = 2( D . (18) 

All exponents cited above are temperature in- 
dependent. Some glassy systems do, however, 
present temperature-dependent exponents asymptoti- 
cally. Whether this behavior can occur in our system 
is a delicate issue, that we discuss in Sect. [V] 

2. Comparison to an aged configuration 

In simple cases in which equilibrium is relatively rapidly 
reached the initial condition should be irrelevant after 
the equilibration time. The same relaxational behavior 
is then expected, independently of the waiting-time t w 
and the initial condition. The roughness should only de- 
pend on the time-difference At = t — t w and the same 
functional forms should characterize growth and satura- 
tion. In an out of equilibrium relaxation relatively soon 
after preparation the growth regime acquires a waiting- 
time dependence and the relations (p~0|) - (|T2]) might be 
generalized. 

For each waiting-time the two-time roughness presents 
two regimes as a function of At: it first grows until 
crossing over at t x (L) to saturation at a At indepen- 
dent value. For high working temperature after a short 
transient the memory of the initial condition disappears, 
and the roughness is well-described by the Family- Vicsek 
scaling. For low working temperatures the waiting-time 
dependence remains. For long lines the crossover time is 
long enough to see aging behavior in a sufficiently long 
time-window, such that we can describe it with a scaling 
form. Before saturation the roughness does not depend 
strongly on the size of the line (for sufficiently long lines) 
and the curves can be scaled as 

(^HM.W 2 ^ <* 2 } • (19) 

This scaling form approaches a stationary regime in 
which (w 2 )(t,t w ) ~ £ 2<: (At) in the limit At > t w if 
(w 2 )(u) ~ u 2 ^ for u 3> 1. For even longer time- 
delays such that £(At) — > L one finds saturation at 



(w 2 )(t,t w ) — > (w 2 ^) ~ L 2 '', which might indeed contain 
a independent prefactor, as in the clean EW caseJ^ 

The out of equilibrium regime exists without the need 
of quenched randomness. We showed in Ref. [H that 
the roughness in the simple EW equation satisfies the 
scaling form ()19|) when the system is let evolve from an 
out of equilibrium initial condition and the waiting-time 
dependence is kept explicit. The length scale £ grows 
as £(t) ~ i 1 / 2 . A similar scaling was shown numerically 
in Refill for the (1 + 1) KPZ model with £{t) ~ t 2 l z . 
In Sect. [TV] we discuss the scaling of Eq. (TT5]) for the 
disordered model. 

The analysis of the saturation of the two-time rough- 
ness should, in principle, yield the thermal and disorder- 
dominated values of £. Computing £ from (w 2 ^) is, how- 
ever, numerically heavy, since one needs to simulate lines 
with different lengths and then extract the scaling be- 
havior. In Ref. [l5| we tried such a scaling analysis but 
we did not reach the regime in which ( = Cd = 2/3 (see 
Fig. 1(d) in this reference). A more convenient way of 
getting C is to study the two-time structure factor^ as 
we explain below. 

B. Structure factor 

1. The roughness exponents 

The thermal, Ct, and disorder, Cd, roughness expo- 
nents can also be extracted from the analysis of the struc- 
ture factor defined in Eq. ([3]). Working with the struc- 
ture factor is convenient since it is sufficient to simulate 
a long chain and then extract the L dependence from 
the wave- vector dependence X Indeed, the roughness is 
simply related to the structure factor as 

oo 

(w 2 ){t } t w ) = L- 1 ( S n)(t,t w ), (20) 

n— — oo 

and this implies that the dynamic scaling jHjJ-fTJ]) (for 
t w — 0) is equivalent to 

(S n )(At) = (S n )(t = At,t w = 0) = S eq (q n ) g[q n £(At)} 

(21) 

with S eq the equilibrium wave-vector dependent struc- 
ture factor and g(u) ~ u d+2< ^ for m«1 and g(u) — > const 
for u 1, with d the internal dimension of the manifold 
(d = 1 for the directed line). These limits for the function 
g(u) ensure, respectively, the equilibrated steady geome- 
try and the memory of the flat initial condition. In the 
case of power-law scaling Eq. (|2Tj) becomes 

(S n )(At) = q- id+2 ° g[q n Ai 1/Z ] . (22) 

In the out of equilibrium relaxation process we are in- 
terested, (S n ) depends on two-times t and t w . The gen- 
eralized scaling form then reads 

(S n )(t,t w ) = S e<1 (q n ) G[q n £(t),q n £(t w )} (23) 
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FIG. 2: (Color online) The value of the wave vector qr signal- 
ing the crossover between the thermal and disorder regimes 
at different working temperatures. The temperature depen- 



dence is consistent with the expected qr 
see Eq. (p4j) . 



T power-law, 



FIG. 1: (Color online) Structure factor in the continuous EW 
line in a random environment for To = 5 and two working 
temperatures, T = 0.5 (circles) and T = 0.8 (squares), with 
t w = 0. The data show the thermal-regime at very large 
wave vector, a crossover to a disorder-dominated regime upon 
decreasing the wave-vector value and, finally, a saturation 
regime demonstrating that the growth- length, £(t), is shorter 
than the system size, L. In the bottom panel the structure 
factor is scaled by q 2 in order to highlight the difference be- 
tween the two roughness exponents. 



that reduces to Eq. (f2Tj) both for t w = 0, and to a similar 
form for q n l{t w ) 3> 1 for all n, and thus stationarity is 
recovered. For fixed t w and At we can study the depen- 
dence of (S n ) on q n . If the aim is to determine the values 
of the roughness exponent the ideal choice would be to 
use a very long At = t — t w ^> t w and plot (S n ) for sev- 
eral t w 's as a function of q n . This construction is to be 
compared with Fig. 2(b) in Ref. 18. In presence of dis- 
order the construction should have a broken straight-line 
form with two exponents, — (1 + 2C,t,d)', the thermal one, 
£t, characterizing the behavior at short length-scales, 
<Zn > qr, arid the disorder one, £d, characterizing the 
behavior at long length-scales, q n < qr- The breaking 
point qr defines a characteristic length Lt ~ 1 /qr which 
should depend on T and we shall discuss and evaluate its 
precise temperature dependence in Sect. 1111 B 21 In an 
infinite (1 + l)-dimensional system disorder always dom- 
inates the fluctuations at very long length scales. How- 
ever, in a finite system of length L once Lt = L all the 
system is characterized by thermal fluctuations and the 
value of Lt is thus crucial. 

Studying the wave- vector dependence of the structure 
factor is thus a convenient way to determine the rough- 
ness exponents, C,t,d] we hence expect to improve the 



values shown in Fig. 1(d) in Ref. [I3. We studied the 
t w = structure factor in the continuous disordered 
EW line. The general behavior described in the previ- 
ous paragraphs is shown in Fig. [T] for working tempera- 
tures T = 0.5 and T = 0.8. Note that the saturation of 
the structure factor at very small wave vector bears the 
same information as the fact that the roughness (w 2 ) has 
not yet saturated at the longest time shown in Fig. [T] 
It means that the growing correlation length has not 
reached the size of the system, £(t) < L. The crossover 
from the thermal, (r, to the disorder, (p, roughness ex- 
ponents is clear, with the expected values Ct = 1/2 and 
(x> = 2/3 confirmed. Note that the top panel in Fig. [1] 
shows the raw data while the bottom panel shows the 
structure factor multiplied by q 2 , which allows to better 
observe the crossover region and the difference between 
the two exponents. 



2. The crossover length scale Lt 

The analysis of the structure factor also allows us to 
evaluate the temperature and disorder strength depen- 
dence of the crossover length Lt (or wave- vector gy) and 
compare it to Yoshino's unpublished results for the lattice 
model^ and Nattermann, Shapir and Vilfan's estimate.^ 6 

The temperature-dependent crossover length Lt sepa- 
rating the saturation regime characterized by the thermal 
and disorder exponents £t and C,d should be observable in 
this presentation as a temperature-dependent crossover 
wave- vector qT in the equilibration prefactor S eq . In or- 
der to access it one then needs to use a sufficiently long 
time-delay in such a way that the second factor saturates 
to g(u) — > const. 
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The crossover length, and the crossover wave-vector, 
are expected to scale with temperature as Lt ~ T 1 ^, 
and qx ~ T -1 /^, with <j> an exponent that we shall ob- 
tain below using a variety of arguments. The charac- 
teristic length-scale Lt ~ should increase with in- 
creasing temperature {<p > 0), indicating that thermal 
fluctuations become more and more important. 

A simple argument to obtain the <p exponent relies 
on the fact that the fluctuations at large length scales 
{in < It) a re n °t expected to strongly depend on the 
working temperature, in contrast with the short length 
scale ones. This is indeed confirmed in our numerical 
simulations, as can be observed comparing the two sets 
of data in Fig. Q] We will use this observation as a 
main input for the scaling argument in Sect. [Vj Let 
us assume that a weak temperature dependence is per- 
mitted, i.e. (S q ) ~ a(T)q n ( 1+2 ^-°) for q n < qT, while 

(S q ) ~ Tqn (1+2CT) for q n > q T . Matching the crossover 
between these two regimes at qT one obtains 

r rp I -1/[2(Cd-Ct)] 



Using a{T) ~ a yields q T ~ T-VPCCu-Cr)] ar fo = 
2(Cd — Ct)- With simulations of the discrete model 
Yoshino found <fi = 1/3/4 in agreement with this result 
since for our one dimensional case <p = 2(Cd — Ct) = 
2(2/3-1/2) = 1/3. 

A different prediction for the exponent <f> comes from 
an order of magnitude argument. If one assumes that the 
characteristic free- energy barrier, A(L), associated to the 
length-scale Lt should be the thermal one, 

~k B T, (25) 

one finds 

L T -L ^y\ (26) 

If T and Lq are independent of temperature, Lt ~ T 1 ^ 
and <f> = tp. Yoshino used this kind of argument with 
the free energy barrier replaced by the sample-to-sample 
fluctuations of the energy, an assumption tested in Ref.l29l 
for low temperatures. In other words he replaced tp by 
9, with 9 the energy exponent to get <p — ip ^ 9 (9 = 1/3 
in d = 1),— This argument assumes that the scale Lt 
is the onset for thermally activated motion, making a 
connection between the dynamics and the geometry of 
the line. 

On the other hand, using Flory-type arguments to pre- 
dict the roughness exponent, Cf = [5 — (d+N)]/5 = 3/5 
for d = 1, with (w 2 ) ~ L 2< > F , and an expansion around 
T = 0, Nattermann et a$&- suggested (p — 2(Cf — Ct) and 
using Ct = [3 - (d + N)}/2 = 1/2 they found 0-1/5 in 
d=N = l. 

Figure [3] displays the result of the analysis of the 
crossover between the thermal and the disorder regimes 
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FIG. 3: (Color online) Logarithmic fit of the roughness re- 
laxation of the elastic string with t w = 0. All the curves 
correspond to the same initial temperature To = 5 and differ- 
ent quench temperatures T = 0.5, 0.7, 0.9, 1.1 and 1.3 from 
lower to upper curves. Circles are raw data and blue contin- 
uous curves are the fitted functions. 

in the t w = structure factor at different working tem- 
peratures. Our numerical simulations of the disordered 
EW line is compatible with <j> — 1/3 in agreement with 
the first two exposed arguments and Ref.[lJ, see Eq. ([24]) . 
but it rules out the value given by Nattermann et al2& 
However, by replacing the Flory approximation to the 
roughness exponent used by these authors by the exact 
exponent, we get 0=1/3 again, in agreenment with our 
results. 



3. The growing length: logarithmic growth from t w = 
measurements 

At high temperature, T > T co , equilibrium is quickly 
reached, disorder is irrelevant and one recovers stationary 
dynamics with the clean EW temperatue-independent 
values a — 1/4, C = Ct = 1/2 and 2 = 2. At low tem- 
perature one has to distinguish between the thermal and 
disorder regimes and analyze the width of the crossover 
region. 

We shall first take t w = and compare with the re- 
sults obtained by Kolton et a@2- (see also Ref. l3l| ) . These 
authors showed, by studying the evolution of the struc- 
ture factor of an elastic line in random media from a flat 
initial condition and using t w = 0, that the characteristic 
growing length crosses over from a power-law dependence 
£(t) — i 1 / 2 typical of the clean EW case to a logarithmic 
law, £(t) - (lnt) 1 /^ with if) - 0.49 (see Fig. 2 in Ref.©. 
The exponent tp characterizes the scaling of the barriers 
at large scales, A(L) ~ L^ . Kolton et al expected ip — 9, 
with 9 the exponent characterizing the free-energy cost of 
an excitation of length L, AF(L) — L e . The value of 9 is 
known exactly for d = N = 1, 6 = 1/3, and the authors 
of Ref. 30 ascribed the discrepancy between the measured 
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FIG. 4: (Color online) Dependence on the working tempera- 
ture of (a) the Qin exponent obtained from the fits in Fig [3] 
and (b) the barrier exponent ip — 2C,o/ot\ rl _, see Eq. (118[) . with 
the dotted lines indicating the suggested values ip = 1/2 and 
ip « = 1/3. The temperature of the initial condition is 
To = 5.0. The error bars were estimated from the dispersion 
of the exponents when changing the fitting range. 



FIG. 5: (Color online) Exponents obtained from a logarithmic 
fit of the roughness relaxation of the elastic string for fixed 
working temperature T — 0.5 and different initial conditions 
Tq. (a) ai n exponent and (b) barrier exponent, ip. The dotted 
line indicates the suggested ip = 1/2 and ip = 1/3 values. The 
error bars were estimated from the dispersion of the exponents 
when changing the fitting range. 



growing length, £(t) ~ (hit) 1 / - 49 , and their expectation, 
£(t) ~ In 3 t to strong logarithmic corrections. 

If the structure factor scales as in Eq. (|2"Tj) , the rough- 
ness scales as 

(w 2 )(At) ~ J d qnq -J 1+ ^g[ qn e(At)} ~ e^(At) . (27) 

Using a long line, L — 256, as in Ref. 30j we confirm 
the logarithmic growth of the roughness for the t w = 
case, as shown in Fig. [3J Different curves correspond 
to increasing temperatures, T = 0.5, 0.6, 0.7, ...1.4 from 
bottom to top. We fit these curves to (w 2 )(At) = 
yl[ln(At/_B)] Qln obtaining the temperature dependent 
a\ n (T) exponent shown in Fig. [Ha). Now, by simply 
using the relation (w 2 )(At) - £ 2<: (At) ~ (InAi)" 1 " one 
recovers the exponent ip, i.e. ip = 2£d / 'a\ n {T), where we 
assumed that in this regime the disorder roughness ex- 
ponent is the relevant one. We obtain ip(T — 0.5) « 0.51 
which is close to the value ip w 0.49 obtained in Ref. H^- 
Recently, Monthus and Garel 3 ^ conjectured ip — d s /2 
with d s the dimension of the surface of the excitation, 
d s = 1 in our case, that is in very good agreement with 
our numerical results at T — 0.5 without any need to ad- 
vocate for logarithmic corrections. However, the analysis 
at different working temperatures presented here reveals 
that the value of ip departs from 1/2 and approaches 1/3 



for increasing temperature, see Fig. H^b). This is the 
contrary to what is expected from the argument given in 
Ref. that is based upon the role of the entropic contri- 
bution. Indeed, one would have expected the value 1/2 
to prevail at high temperatures, which is the opposite 
trend to what we find in the numerical simulations. We 
shall show in Sect. [V] that the temperature dependence 
of ip observed in Fig. HJb) can be explained as a finite 
size/time effect induced by the crossover. 



We have also analyzed the behavior of these exponents 
when changing the initial temperature, while keeping 
t w = 0, that is to say, using initial conditions thermalized 
at different Tq's. We performed the same fitting proce- 
dure keeping the working temperature fixed to T = 0.5 
and varying T = 2.0, 2.5, 3.0, 3.5, 4.0, 4.5, 5.0. The re- 
sulting exponents are shown in Fig. [SI Despite the rather 
large fluctuations we can trust that there is no system- 
atic dependence on the initial temperature. We conclude 
that the exponent values do not depend on the previous 
history but they do on the relaxation conditions. 
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4- The growing length: pre- asymptotic power-law at finite 



When the waiting time dependence is considered, the 
logarithmic growth is no longer easily observed. On the 
one hand, the crossover between the two asymptotic lim- 
its, t w <C At and t w 3> At, makes it difficult to fit a 
logarithmic growth. On the other hand, when the corre- 
lation length £ increases the necessary time to reach the 
logarithmic regime increases exponentially r ~ exp(^). 
Thus, from a practical point of view, the logarithmic 
regime is not reached and one can use a power-law growth 
of the roughness as an effective description. We shall use 
such a power-law to analyze the aging behavior of the 
roughness in Sect. IIV Al 

In this case, the best option to extract the dynamic 
growing length £(t) is to take a very long t w and study 
the breaking point q x ~ \ jt in the structure factor. This 
breaking point separates the modes keeping the memory 
of the initial condition, q <C q x , from the equilibrated 
modes described by the power-law regime with disorder 
exponent, q q x . In the case of a power-law scaling, 
t ~ L z , one has 

q x ~ Ar 1 / 2 (28) 

and from here one obtains z. Using power-law fits we 
find that the growth and dynamic exponents, a and z, 
acquire a T-dependence in the presence of quenched dis- 
orde r 12 ' 13 ' 14 ! 15 while the roughness exponents Ct.d do not 
depend on temperature. 



IV. AVERAGED AGING DYNAMICS 

In this Section we analyze the two-time evolution of 
quantities that are averaged over many realizations of 
the thermal noise and quenched disorder. We consider 
two types of preparation: cooling from a higher temper- 
ature and heating from zero temperature. We discuss 
how the initial condition affects the behavior in the pre- 
asymptotic growth regime. 



A. The two-time roughness 

The analysis of the averaged two-times roughness of 
a very long line described by the solid-on-solid model 
followin g a qu ench from infinite temperature was given 
in Refs. lT2lll3lfl4lll5l . We summarize here these results 
and we complement them by showing that: i. a simi- 
lar scaling form describes the relaxation dynamics of the 
disordered EW line; ii. the dynamics after heating up 
the lines can also be described by the same scaling form 
with parameters that depend on the working and initial 
temperature; Hi. the qualitative behavior is similar to 
the one found analytically for the clean E W 18 ' 18 and nu- 
merically for the clean KPZ22, lines. 
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FIG. 6: (Color online) Panel (a): the noise and disorder aver- 
aged two-time roughness of the lattice elastic line with length 
L = 500 evolving atT = 5, 0.5, 0.2 after a quench from in- 
finite temperature, To — * oo. Panel (b): scaling plot of the 
averaged two-time roughness of the line with L = 500 and 
T = 0.2, and the line with L = 5000 and T = 0.5. 



The time-difference dependence of the averaged rough- 
ness evolving from an infinite temperature initial condi- 
tion in the SOS model, for several waiting-times and at 
three working temperatures, is shown in Fig. Ela). The 
averaged roughness of the line evolving at high temper- 
ature, T — 5, is stationary. It grows as the clean EW 
roughness, (w 2 ) ~ At 1 / 2 and it reaches saturation at 
(w 2 ^) after At > t x . The line with L — 500 evolving at 
T = 0.5 is not stationary but signs of saturation are vis- 
ible at At ~ 10 6 MCs. The line evolving at low temper- 
ature, T = 0.2, shows aging effects, is still in the growth 
regime and does not reach saturation in the numerical 
time-window. 

The averaged roughness curves in the growth regime 
can be well described by the scaling in Eq. I|19p. The 
numerical analysis of the growing length £(t) in the time- 
window we can reach indicates that £(t) ~ t x l z . As 
discussed in Sect. IIIIB 31 this power-law regime should 
be considered as a pre-asymptotic approximation to a 
slower, logarithmic growth. We shall continue our anal- 
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FIG. 7: (Color online) Temperature dependence of the effec- 
tive growth exponent a and the parameter B measuring the 
distance between the two asymptotes Co and c x in the lattice 
model evolving from an infinite temperature initial condition. 



ysis restricting to the power-law regime. 

The resulting f/^-dependent factor (w 2 ), see Eq. (|19p . 
can be well acquainted for by the empiric formic 



r.2N ( 



g(x) 



with 



x — At/t v 



G(X) ~ X 2a ^ A(T, T ) 1Q B(T,T 0)g ( x ,T,T ) (2Q) 

( X 



g(x,T,T ) = tanh 



C(T, T ) lo gl 



310 \D{T,T ) 
Equations (fT9|) and ((|29|) have the stationary limits 



(w 2 ) 



c Q (T,T ) At 2a m At«t„ , 
Coo (r,r ) At 2 "( T ) At»t U(> 



(30) 



with 



c (T,T ) = A(T,T ) 10~ s ™ 
Coo (T,T ) = A(T,T ) 10 S ^ T °) . 



(31) 



These parameters control the two stationary asymptotes 
bounding the aging regime. The parameter 2B(T,T ) = 
log 10 [coo (T, Tq)/co(T, T )] is then a measure of the dis- 
tance between the two asymptotes. 

Figure [S^b) shows the scaling plot of the averaged 
roughness in the growth regime. We display data for 
L = 5000 at T = 0.5, and L — 500 at T — 0.2. Note 
that the data for L — 500 at T = 0.5 shown in panel (a) 
of Fig. [S] would have not scaled properly since saturation 
appears earlier. 

The temperature dependence of the (effective) growth 
exponent a as well as the parameter B for the solid-on- 
solid model evolving from an infinite temperature ini- 
tial condition are shown in Fig. [71 The exponent a ap- 
proaches the EW value 2a = 1/2 at high temperature. 
In the high-temperature limit the parameter i?(T, oo) 
approaches zero since the two constants cq(T, oo) and 
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FIG. 8: (Color online) The averaged two-time roughness in 
the continuous model after a quench from high temperature 
To = 5. (a) L — 256 and different temperatures. The up- 
per curve corresponds to the evolution of the roughness when 
the temperature is the same as the initial value, i.e. a sta- 
tionary situation. Curves for T = 1 (middle curves) and 
T — 0.5 (lower curves) are shonw. The waiting times are 
t w = 100,300,1000 and 3000 from left to right, (b) Scaling 
plot of the data presented for T = 1.0 and T = 0.5 using the 
a values given in the figure. 



C 




FIG. 9: Initial and working temperature dependence of the 
scaling exponent a in the continuous EW model in a random 
environment, (a) The initial temperature is fixed to To = 5 
and the working temperature is varied, (b) The exponent a 
does not change while changing the initial temperature To and 
keeping the working temperature fixed at T = 0.5. 
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Coo(T, oo ) tend to the same temperature independent 
value and aging disappears. 

Figure [5] shows the averaged two-time roughness for 
the continuous model with a high initial temperature 
T = 5. The upper curve in Fig-JSJa) corresponds to the 
stationary case where the temperature is kept constant 
at the initial value. The other two cases, with T = 1 and 
T = 0.5, present aging which is qualitatively similar to 
the SOS model of Fig. [6] FigureOJb) shows the scaling of 
these curves with the temperature dependent exponents 
a given in the key. In the continuous model the expo- 
nent a increases with temperature until reaching the EW 
value 2a = 1/2 after the crossover temperature T co . In 
Fig. [Ha) we present a(T) for the cases with a fixed initial 
temperature To = 5. As can be observed, its precise tem- 
perature dependence is not easy to determine. Since one 
expects that the a exponent determines the slow dynam- 
ics at a given temperature, it is interesting to check that 
the value of a at a small temperature does not depend 
on the initial temperature. This is shown in Fig. [9jb), 
where we show that a(T = 0.5) is independent of the 
initial temperature. Although we do not compute the 
parameter B for the continuous model, it is worth noting 
that in contrast with the discrete case the two constants 
Co and Coo approach a temperature dependent asymptote 
Co ~ Coo oc 2T that equals the high temperature limit of 
the EW line, c oc 2T and cx (2- V2)T + V2T . 18 The 
exponent a monotonically increases from a(T = 0) ~ 
and crosses over to 2a — 1/2 at some T co with a very 
weak T dependence. 

The effect of using as an initial configuration one in 
equilibrium at a lower temperature was considered ana- 
lytically in Ref. [l8| for the EW line. In presence of disor- 
der one is often forced to use numerical simulations and 
the difficulty of equilibrating a long line at low tempera- 
ture arises. However, a special case is that of the lattice 
model at To = 0. It is indeed well known that its ground 
state configuration can be exactly calculated by transfer- 
matrix methods j2£ 

In Fig. [TU] we show the two-time roughness of the dis- 
crete disordered line using as initial configurations equi- 
librium ones at To = and a working temperature 
T = 0.2. The plots demonstrate that the trend of the 
curves is to reverse in the sense that the curves go from 
the upper to the lower asymptote. This behavior is simi- 
lar to what was found analytically for the clean EW line 18 
and the KPZ non-linear model^S This effect is also rem- 
iniscent of what was observed in the relaxation of the 2d 
XY model from a uniform initial condition, correspond- 
ing to equilibrium at To = 0, at finite (higher) tempera- 
ture within the spin-wave approximation^ 

The effect of the initial condition can be summarized 
in the following way. At all working-temperatures that 
are identical to the one of the initial condition, T = To, 
one finds Co = Coo and there is no aging since the line 
is initially in equilibrium^ For T > T one has c > 
Cqo while for T < T one has Co < Coo- The values of 
these constants determine the relative location of the two 
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FIG. 10: (Color online) Panel (a): the averaged two-time 
roughness in the lattice model for an elastic line with length 
L — 5000. Evolution at T = 0.2 after a sudden heating from 
equilibrium at To = 0. The waiting times are tw = 10 (plus, 
red), t w — 10 2 (cross, green) t w = 10 3 (star, blue), t w — 10 4 
(open square, pink), t w = 10 (filled square, cyan). Panel (b): 
scaling plot of the data in panel (a). Note that the value of 
a (in the key) coincides with the one used in the case of a 
quench to T from To —> oo (see Fig. 0. The dashed lines are 
the bounding asymptotes. Let us remark that, here, Coo < Co- 

asymptotes of the growth regime, At -C t w and At ^> t w . 

B. The linear response 

We now turn to the study of the linear response. In 
Fig. [TT] we show the averaged integrated linear response 
of the roughness in the lattice model. In panel (a) we 
show data for several waiting-times as a function of time 
delay. In panel (b) we scale the data as 

(*>(*,*«) ~^(*«)<X>(^) (32) 

with l{t) ~ t 1 ' 2 . 

These results are to be confronted with the behavior 
of the EW elastic line with Langevin dynamics for which 
the integrated linear response is stationary and does not 
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FIG. 11: (Color online) The noise and disordered averaged 
integrated linear response associated to the two-time rough- 
ness of the lattice elastic line. Panel (a): the line with length 
L = 500 evolving at T = 5, 0.5, 0.2 after a quench from in- 
finite temperature, To — *■ oo. Panel (b): scaling plot of the 
averaged two-time linear response of the line with L — 500 
and T = 0.2, and the line with L = 5000 and T = 0.5. 



fluctuated On the other hand, the linear response ob- 
tained with the continuous model, not shown here, is 
qualitatively similar to the one shown for the SOS model 
in Fig. HU 

In Fig. [12] we display data for the linear response af- 
ter the sudden heating from the ground state used in 
Sect. lIV A1 The effect of this 'reversed' heating-procedure 
is similar to the one showed in Sect. IIV Al and in Fig. [TO] 
for the roughness. 



C. The Fluctuation-dissipation theorem (FDT) 

The modification of the FDT linking the displacement 
to its associated response in the lattice model after a 
quench from infinite temperature was studied by A. Bar- 
raid and Yoshino i 13 ' We here focus on the behavior of 
the roughness and its linear response. We show data for 
the disordered EW line after a similar quench. In Fig. [TBI 
we show the plot (x) against (w 2 ) at fixed t w and using 



FIG. 12: (Color online) Panel (a): the averaged integrated 
linear response associated to the two-time roughness in the 
lattice model for an elastic line with L = 5000. Evolution at 
T — 0.2 after a sudden heating from equilibrium at To = 0. 
The waiting times and symbols are as in Fig. 1101 Panel (b): 
scaling plot of the data in (a). Note that the value of a (in the 
key) coincides with the one found in the study of (w 2 ) and the 
one used in the case of a quench to T from To — * oo (see Figs. [6] 
and llO|) . The dashed lines are the bounding asymptotes. Let 
us remark that, here, Coo < cq. 



At as a parameter going from At — to At — > oo in the 
continuous disordered EW line for the case Tq = 5 and 
T = 0.5. As it can be observed, it displays two slopes, 
allowing for the definition of an effective temperature 7 . 
A similar behaviour was reported for the displacement 
and associated linear response in the lattice modeld the 
clean EW lined and the vortex glass models 

The temperature dependence of the effective temper- 
ature is shown in Fig. [U for the disorder and the clean 
case. The dependence of T e g on the working temper- 
ature while keeping To fixed is shown in Fig. [JJJa), 
but it is difficult to asses whether the effective tem- 
perature is constant or slowly grows with T. With- 
out disorder the EW solution gives a linear dependence 
T eS = T/V2 + T {1 + 1/V2), 18 shown with dashed lines 
in the figure. When the working temperature is fixed 
at T = 0.5 and the initial temperature is changed, T e g 
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FIG. 13: (Color online) The parametric plot of the scaled 
linear response against the roughness for the continuous model 
with To = 5 and T = 0.5. The dashed line indicates the FDT 
limit. 
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FIG. 14: (Color online) Temperature dependence of the effec- 
tive temperature, signaling the violation of the FDT, compar- 
ing the clean EW system -dashed lines- with the disordered 
continuous case, (a) Dependence on the working temperature 
for a fixed initial temperature To = 5. (b) Dependence on the 
initial temperature To while the working temperature is kept 
fixed at T = 0.5. 



grows linearly with To, as shown in Fig. I14f b) and in the 
clean limit. 

We also studied the effect of a heating procedure on the 
violations of the FDT using the lattice model. We show 
in Fig. [T5] the parametric plot (x) against (w 2 ) using as 
initial condition the ground state at T = 0. As already 
found in the clean EW line we find that T c g < T. T c ff 
thus reflects, consistently, a 'memory' of the initial con- 
figuration. 



V. CROSSOVER-INDUCED GEOMETRICAL 
AND DYNAMICAL EFFECTS 

We discuss here how our numerical results can be ex- 
plained by developing simple scaling arguments based on 
the existence of a single dynamic crossover in the grow- 
ing correlation length at a static temperature dependent 
length Lt, from a thermally dominated regime, to a dis- 
order dominated regime with algebraically growing bar- 



FIG. 15: (Color online) The parametric plot 2T(x) against 
(w 2 ) after a heating procedure from the ground state at To = 
to T = 0.2 in the SOS model. The waiting times are 10 MCs 
(plus, red), 10 2 MCs (cross, green), 10 3 MCs (star, blue), 10 4 
MCs (open square, pink) and 10 s MCs (filled square, cyan). 
Note that T eB <T. 



riers as function of the length-scale. The main additional 
assumption on that we make is that at large enough 
length-scales I Lt, the geometry and typical barriers 
controlling the dynamics are indistinguishable from those 
at zero temperature. With these hypothesis, which are 
qualitatively supported by our numerical results, we can 
predict: (i) the temperature dependence of Lt; (ii) the 
temperature dependence of the effective exponents char- 
acterizing each regime of growth; (iii) the temperature 
dependence of the crossover time between these regimes; 
(iv) a parameter quantifying the importance of crossover- 
induced finite-size/time effects in simulations. In order 
to make the discussion in this Section self-contained we 
start by repeating some arguments that are well-known 
in the literature and then we present the scenario that, 
we propose, explains the numerical observations. 

If the temperature is high enough to renormalize the 
microscopic pinning parameters^ thermal fluctuations 
dominate the line wandering at short length-scales. The 
geometry of equilibrated small scales, I <C Lt, is thus 
described by the thermal roughness exponent (t = 1/2 of 
the clean Edwards- Wilkinson (EW) equation, such that 
(w 2 ) ~ (T/c)£ for I/Lt < 1 or S(q) ~ (T/c)g-( 1+a ^) for 
qLt S> 1, or q 3> <7t ~ The typical time needed to 

equilibrate a length £ in this regime is therefore expected 
to be t(£) ~ £ z , with the EW dynamical exponent z = 2. 

At length-scales £ 3> Lt, disorder dominates and the 
geometry of the equilibrated large length scales, £ 3> Lt, 



is the same as in the ground-state, S(q) 



^ -U+2Cd) 

for qLT <§C 1, with £d = 2/3 > Ct- The typical time 
to equilibrate a length I ^> Lt is controlled by size- 
dependent energy barriers U(£), such that t(£) ~ e u ^^ T , 
with U{£) ~ ^> T. By hypothesis, we assume that 
both the form of S(q) and the barriers U{£) are indepen- 
dent of T in this regime, in agreement with what is seen 
in Fig. [T] for the structure factor with q <C gr- 
in order to capture the crossover effects we interpolate 
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the geometric and dynamic behavior of the regimes I 3> 
Lt and I -C Lt described above. The static structure 
factor can be written as 



S(q,T)~q-^+^g(qL T ,T) 



(33) 



with g{x,T) ~ Tx 2 ^°-^ for x > 1 and g(x,T) ~ 
const, for This interpolating form, which includes 

the temperature, is found to fit well the data, and to 
satisfy our assumption S(q, T) ss S(q,T = 0) for qLr <C 
1. We continuously match the two regimes by requiring 



Tq, 



-(l+2Cr) 



and, therefore, 



(34) 



with <p = 2(£d — Ct) = 1/3, which is the temperature 
dependence of the crossover length we observe in Fig. 
[2j This argument was also discussed in Sect. MI B 21 on 
numerical grounds, where references were also given. 

Let us now discuss the dynamical behavior induced 
by Lt- At the crossover we can roughly set U(Lt) ~ 
T, since barriers start to dominate the dynamics only 
above L T . For £ > L T we expect U(£) ~ & > T. To 
interpolate the large scale behavior down to Lt we use=£ 



U(£) 



W/Lt)*. 



(35) 



Since by hypothesis U{£, T) w [/(.£, T = 0) for I > L T , 
we obtain ip — <p — 2(Cd — Ct) = 9 = 1/3, with the 
energy exponent. The time to equilibrate a length £ can 
hence be written as 



t{£) 



(36) 



This expression continuously matches the barrier- 
dominated regime, t ~ exp[(£/ Lt)^] = exp[£^/T] for 
I 3> £tj with the power-law growth expected in the ther- 
mal regime, t ^ l z for £ <C L^. From Eq. (|3"6"|) and using 
Lt - T 1 /^ one obtains a characteristic time scale, 



T 



z/lp rpQ 



(37) 



separating thermally from disorder dominated dynam- 
ics. This strong temperature dependence explains why 
the crossover to the slow logarithmic growth is partic- 
ularly difficult to observe in numerical simulations. On 
the one hand, at high temperature, the line size must 
be large L > Lt ~ T 3 in order to see the crossover to 
the long-time logarithmic regime. On the other hand, at 
low temperature, the lines could in principle cross over to 
the logarithmic regime even for small system sizes but, 
because of the very slow Arrhenius activated dynamics, 
exceedingly long running times would be needed to re- 
solve the precise time dependence. Most importantly, the 
data from any numerical simulation displaying the two 
regimes of growth will be characterized by temperature- 
dependent effective exponents due to crossover-induced 
finite-size/time effects, as we explain below. 

Indeed, Eq. (|36[) allows us to explain the temperature 
dependence in the effective exponents of the power-law 



and logarithmic growth regimes. The influence of the 
crossover in the effective dynamical exponent z of the 
power-law regime relevant to short length-scales, £ -C Lt, 
follows from equating Eq. (f3"6")) to t(£) ~ £ z , which yields 



L 



(38) 



To avoid solving the resulting self-consistent equation, we 
can use z ~ z in the third member to obtain a first order 
correction to the dynamical exponent, 



(2^/7 In i) 



(39) 



Since by construction the power-law fit of the data is 
done in a relatively short time-scale such that t^' z / hit is 
almost constant (otherwise z would not be well defined), 
and by using Lt ~ T 3 and ip = 1/3, we get 



(40) 



This rough estimate is in good agreement with the nu- 
merically determined exponent a shown in Figs. [7] and [5] 

It is worth noting here that the Cardy-Ostlund ran- 
dom Sine-Gordon model 34 in d = 2, disordered ferromag- 
nets,— and Ising and vector spin glasses^ were shown to 
display a dynamic exponent z ~ 2 + a/T, with a of the 
order of unity and depending on details of the considered 
model. Our arguments and numerical study show that 
the numerically measured dynamic exponent could sim- 
ply be an effective value that stems from the interpolation 
of two growth regimes that cannot be properly resolved 
numerically. This seems to be particularly important for 
systems with a positive energy exponent 9, since in gen- 
eral we expect ip ~ 9, and therefore a crossover from 
a thermally dominated short-time regime to a long-time 
dynamics dominated by diverging energetic barriers. 

The temperature dependence of the effective exponent 
in the regime of logarithmic growth, £ 3> Lt, can be 
explained with a similar argument. By equating Eq. (|36[) 
to t{£) ~ exp[f^/T] we find, 



ip — ip 



\n[T/L T + zT In lj£^\ 



hx£ 



(41) 



We use now Lt ~ T 3 and the time-dependence £ ~ 
(Tlni) 1 / 1 ^ in the second member to get a self-consistent 
equation for ip, namely 



ip — ip 



iP 



In (Tint) 



In 



Tz ln(Tlnt) 



(42) 



Fixing the time scale and setting ip sa ip in the second 
member we obtain a first order approximation for ip — ip, 



ip — ip 



iP 



ln(Tlni) 



hi 



zln(Tln<) 
ip hit 



(43) 
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At high temperatures (or large times) such that T In t 
1 this expression simply yields %p ~ ip. In the opposite 
low T limit for sufficiently large length scales such that 
the restriction I ^> Lt ~ T 3 is still satisfied one finds 
ip > ip. Indeed, ip is a decreasing function of T. These 
results are consistent with the result shown in Fig. 2J 
where we observe that ip > ip at low temperatures, and 
ip — ► = 1/3 increasing the temperature. 

The effective exponent analysis made above reveals 
that crossover induced effects are controlled by the or- 
der of magnitude of i^/ln£. Since I in a simulation is 
always a fraction of the system size L, a temperature de- 
pendence of the effective exponents can be expected for 
systems such that L 1 ' 3 /lnL ~ 0(1). For instance, in 
order to properly resolve the dynamic exponent z = 2 
at short length scales we need I <C Lt for at least a 
few orders of magnitude range of £, i.e. a large Lt- On 
the other hand, to properly resolve the exponent ip we 
need f-^/hxi > 1 for L > i > Lt for at least a few 
orders of magnitude of i. Typically, the largest systems 
analyzed numerically in the literature have L ~ O(10 3 ), 
with L 1 / 3 / InL ~ 0(1) and do not satisfy the latter con- 
straint. This clearly prompts for a careful interpretation 
of the anomalous temperature dependence of power-law 
or logarithmic growth exponents obtained numerically so 
far, as we discuss below. 

In Ref. the relaxation of an elastic line in a dis- 
ordered potential was analyzed for the particular case 
t w = and T = 0.5. The crossover to the logarithmic 
growth was observed, though with an effective exponent 
ip w 0.49 > 9. In Ref. [32| scaling arguments for the 
free energy of droplet excitations were given favoring the 
value ip = 1/2, against the commonly accepted value 
ip = 9 = 1/3, 28 in apparent agreement with the results of 
Ref. In the present work we find ip ~ 0.51 for t w = 
and T = 0.5, consistent with the result of Ref. H3, but 
also a decrease towards a value close to ip = 1/3 with 
temperature. This behavior is opposite to what is pro- 
posed in Ref. HH, where entropic effects leading to the 
value = 1/2 are expected to grow with T, from the 
T = value ip = 9 = 1/3, to the entropically- driven 
value ip = 1/2. The scaling arguments given above thus 
show that the anomalous temperature dependence of ip 
can still be interpreted without dropping the ip = 6 = 1/3 
identification, which was proposed time ago^ 

In numerical simulations of the steady-state slow 
driven motion of an elastic line— the effective creep expo- 
nent \x was found to be temperature dependent. Accord- 
ing to the assumption U(£) ~ £ e at large length-scales, 
simple scaling arguments lead to fj, — 9/(2 — (r>) = 1/4 
(see, for instance Ref. l43h . Taking into account that 
L^/hxL ~ O(l) in those simulations we can expect, 
following an equivalent line of reasoning as above, a 
temperature dependent effective creep exponent jl = 
■0/(2 — Cd) > such that jj, —> fj, when T grows. This 
behavior is i n g ood qualitative agreement with the re- 
sults in Ref. 1371 Therefore, the anomalous temperature 
dependence observed in Ref. H3 can not be used as ev- 



idence against obtaining the creep exponent entirely as 
a function of static exponents, although the large scale 
geometry in the creep regime is found to be described by 
depinning exponents i 37 ! 38 



VI. CONCLUSIONS 

Elastic manifolds are objects that appear in a large 
variety of problems with glassy features. It was 
shown in the past that these systems present aspects 
of glassy dynamics combined with diffusion proper- 

^ c „ 6 .11.12.13.14.15.18.19.20 

After a series of analytic and numerical studies of con- 
tinuous and discrete models with and without quenched 
disorder we can summarize the main features of the out 
of equilibrium relaxation of individual and interacting di- 
rected one-dimensional objects. 

All these systems age below a characteristic tempera- 
ture T co . In the mean-field limit of an infinite number 
of transverse dimensions, for a line with infinite length, 
a dynamic phase transition arise at T co . Below T co the 
lines age without diffusion 

The dynamics of finite length lines moving in finite di- 
mensional spaces cross over from diffusive- aging growth 
to a regime in which the roughness saturates. This phe- 
nomenon can be described with a generalization of the 
Family- Vicsek scalingi&i^ The qualitative features of the 
two-time freely relaxing observables are generic and do 
not depend on the presence of quenched randomness but 
the details such as the exponents and growing length do. 
We performed a careful analysis of the time-dependent 
growing length and we found a crossover from an effective 
power-law with temperature dependent exponent to the 
expected asymptotic logarithmic growth for sufficiently 
long strings. Our numerical data indicates a temperature 
dependent effective barrier exponent ip, which is higher 
than the energy exponent 9 = 1/3 at low temperatures, 
but tends to it at high temperatures, provided the sys- 
tem is large enough to display the asymptotic behavior. 
As we discussed via scaling arguments this behavior can 
be recast in terms of crossover induced effects, and it is 
not inconsistent with the usual assumption -0 = 9 for this 
system. 

We also found a single crossover length Lt ~ T 3 sep- 
arating a thermally dominated regime with £t = 1/2 
and a disorder dominated regime with the T = rough- 
ness exponent — 2/3. We showed that this is consis- 
tent with the assumption that large scale properties of 
the system are indistinguishable from those at T = 0. 
Translating this crossover into the dynamics we could 
also describe the temperature dependence of the effective 
dynamical exponents in the power-law growth regime, 
and in the logarithmic growth regime. In this respect it 
would be important to study the corresponding crossover 
in (1+2) dimensions, relevant for vortices in superconduc- 
tors. This is crucial for high-Tc superconductors where 
thermal fluctuations are known to strongly renormalize 
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the disorder and produce Larkin lengths growing expo- 
nentially fast with temperature^ As the Larkin length 
marks the onset of metastability and the crossover to a 
barrier dominated (random-manifold) regime, as it does 
Lt in our case, we can expect the phenomenology and 
crossover induced effects we describe here to be experi- 
mentally relevant, specially at the onset of irreversibility, 
near the transition to the vortex liquid phasej^ where 
the Larkin length can become comparable with the sam- 
ple size4i 

The importance of measuring linear responses was 
demonstrated by the fact that all known coarsening 
systems (above their lower critical dimension) have an 
asymptotically vanishing linear response in the aging 
regime, in contrast to solvable mean-field glassy models 
and numerical simulations of finite dimensional glasses 
that yield a finite integrated linear response in the same 
two-times regime (see e.g. Refs. HI and l39h. This fact 
appears as a concrete difference between the relaxation 
dynamics of coarsening and glassy systems. 

The relaxation of the integrated thermal averaged lin- 
ear response of elastic manifolds strongly depends on 
the presence or absence of quenched disorder. Mean- 
field elastic manifold models in quenched random envi- 
ronments have a fully aging linear response. The clean 
EW line has a stationary linear response while the dirty 
continuous or lattice models have integrated linear re- 
sponses with aging and diffusion combined. 

The effective temperature 7 defined from the compar- 
ison of induced and spontaneous averaged fluctuations 
is finite in all the cases considered as long as the diffu- 
sive factors in the correlations and linear responses are 
divided away. We also showed that the effective tem- 
perature is higher than the bath temperature for cooling 
procedures and, inversely, it is lower than the bath tem- 



perature for heating procedures. This is similar to what 
has been previously found in the 2d XY model 33 and 
the Edwards- Wilkinson elastic line^ and gives support 
to the notion of effectvie temperature as measured from 
deviations from the fluctuation-dissipation relation. 

The study of dynamic fluctuations in these systems has 
not been fully developed yet. In Ref. [H we analyzed the 
fluctuations of the two-time roughness of the lattice disor- 
dered model during growth. The equilibrium^ and out of 
equilibrium 18 roughness fluctuations in the EW line show 
a similar pattern. The probability distribution functions 
in all these models satisfy a scaling law and the scaling 
function follows the same trend as a function of all its 
variables. The question then arises as to whether the fluc- 
tuations of the linear responses of clean and disordered 
elastic lines are similar or different. We shall discuss this 
problem in a separate publication. The relation with the 
theory of fluctuations based on time-reparametrization 
invariance— will also be discussed elsewhere. 
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